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We present detailed calculations of the electron thermalization by synchrotron self-absorption accounting for cooling 
by Compton scattering. For the first time, we solve coupled kinetic equations for electrons and photons without 
any approximations on the relevant cross-sections and compute self-consistently the resulting electron and photon 
distributions. The results presented in this contribution may be applied to the magnetized coronae around Galactic 
black hole accretion disks. 
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1. Introduction 

We have developed a numerical code to study ra- 
diative processes in relativistic magnetized plasmas. 
Previous works made various simplifications mak- 
ing them useful only in certain parameter regimes. 
Usually one assumes that Compton scattering takes 
place in the Thomson regime. Often one treats 
all electrons as ultrarelativistic and neglects dif- 
fusive thermalizing processes such as synchrotron 
self- absorption operating at lower energies. In some 
works the kinetic equations are not solved self- 
consistently (see e.g. [1]), instead some assumptions 
are invoked for the distribution of one or the other 
type of particles. One of the well-known codes deal- 
ing with radiative processes in similar conditions is 
eqpair by Coppi [2], which however does not con- 
sider synchrotron thermalization, but instead treats 
Coulomb collisions as the thermalizing process. 

The only major simplifying assumption in our 
treatment concerns the geometry of the active region. 
Our calculations are done in a simple one-zone geom- 
etry with a tangled magnetic field and isotropic par- 
ticle distributions. However, we have tried to avoid 
the other approximations mentioned above by ac- 
curately treating the microprocesses involved. We 
have used the exact Klein-Nishina Compton scat- 
tering cross-sections [3] for both types of particles 
and exact cyclo-synchrotron emissivities. The cou- 
pled time-dependent kinetic equations for electrons 
and photons are solved self-consistently, making no a 
priori assumptions about the distributions of cither 
type of particles. Such approach makes our code ap- 
plicable in a wide range of parameters. 



2. Code description 

We solve numerically coupled kinetic equations for 
photons and electrons: 
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Here the quantities Df /Dt describe the rates of dif- 
ferent processes (C - Compton scattering, S - syn- 
chrotron) and contain both integral and differential 
terms. The quantities Q p h and Q c denote sources of 
particles not associated with synchrotron or Comp- 
ton processes, such as injection of high-energy elec- 
trons and external sources (e.g. disk) of photons. 

The treatment of Compton scattering over a 
wide energy range is not straightforward due to the 
very different behaviour of the process in different 
regimes. We have therefore adopted an approach, 
where we treat the process as discrete or alternatively 
as a continuous energy gain/loss mechanism depend- 
ing on whether or not the redistribution function can 
be resolved by our grid. This results in the presence 
of both integral and up to second order differential 
terms (to account for diffusion/thermalization) in 
Df/Dtc for both photons and electrons. To account 
for synchrotron cooling and thermalization, the rate 
Df e /Dts also contains up to second order deriva- 
tive terms. The resulting integro-differential equa- 
tions (1) and (2) are discretized on energy and time 
grids and solved iteratively as two coupled systems 
of linear algebraic equations. 
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Fig. 1. Equilibrium electron (left) and photon (right) distributions for different injection slopes: r e = 2 (solid line), 3 (dotted 
line) and 4 (dashed line). Other parameters: tt = 2, L = 10 37 erg s _1 , »?b = 5, no external radiation (/ = 0). The Maxwellian 
parts of the distributions correspond to electron temperatures of 12, 24 and 36 keV, respectively. 



3. Setup 

To demonstrate some of the capabilities of our code 
we have arranged a simple setup, where we inject a 
power-law distribution of electrons into the emission 
region permeated by the magnetic field and let them 
cool and thermalize by synchrotron and Compton 
processes. The seed photons for Compton upscatter- 
ing are produced by synchrotron emission, in the last 
figures we also add an external source of soft disk 
photons. The main parameters are as follows: 

a) r o : injection slope of the power-law, 7 _r °, elec- 
trons, extending from 7 = 1 to 100 in all cases. 

b) tt: equilibrium Thomson optical thickness. 

c) R: size of the emission region, fixed at 3 x 10 7 cm. 

d) L: total power entering the active region as elec- 
tron injection and external soft radiation. 

e) /: ratio of external soft photon compactness (di- 
mensionless ratio L/R) and injected electron com- 
pactness, ZdiskAnj- 

f) r/B- ratio between magnetic and electron injection 
compactnesses. 

g) Tbb- disk blackbody temperature (not a free pa- 
rameter, but depends on L and R). 



4. Numerical results 

We now present the results of simulations for differ- 
ent parameters as equilibrium electron and photon 
distributions. 



4.1. Variable injection slope 

The low-energy electrons absorb synchrotron radi- 
ation produced by the high-energy tail of the elec- 
tron distribution. The more power is in the non- 
thermal tail (the flatter the tail) the higher is the 
self-absorption frequency and the more seed photons 
there are for Compton upscattering, in accordance 
with [5] . Note that a larger amount of seed photons 
results in a drop of the electron temperature (Fig. 1, 
left). Also take note of the curious secondary "bump" 
that develops in case of relatively flat injection spec- 
trum. It is caused by the fact that synchrotron emis- 
sion produced by these electrons is still strongly self- 
absorbed, while the energy losses and gains stay close 
to each other for an extended energy interval. In fact, 
from simple arguments one can show that the ratio of 
energy loss and gain rates for relativistic power-law 
electrons emitting in self-absorbed regime is approx- 
imately 



7c 
7h 



P 



(3) 



where p is the index of the steady-state power-law 
electron distribution. Observe that for p = 3 the 
heating and cooling rates are balanced. However, as 
discussed by Rees [4] already in 1967 and as can be 
seen from (3), such quasi-equilibrium is unstable. 

As we see from Fig. 1, the Comptonized spec- 
trum for hard injection T c = 2 is much softer than 
the spectra observed in the hard state of Galactic 
black holes (GBH), even without any contribution 
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Fig. 2. Equilibrium electron (left) and photon (right) distributions for varying injection power: L = 10 36 (solid line), 10 37 
(dotted line), 10 38 erg s — 1 (dashed line). The electron escape timescale is held fixed, resulting in different optical thicknesses in 
the steady-state: tt = 0.2 (solid), 2 (dotted), 20 (dashed), respectively, with corresponding equilibrium electron temperatures of 
140, 30 and 1.3 keV. Other parameters: Y e = 3.5, rfs = 5, no external radiation. 
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Fig. 3. Equilibrium electron (left) and photon (right) distributions for variable Vb- Vb = ^ (solid), 5 (dotted), 10 (dashed). Other 
parameters: r e = 3.5, tt = 2, L = 10 37 erg s -1 , no external radiation. Equilibrium electron temperatures are around 30 keV. 



to the cooling by the disk. This strongly constrains 
the electron injection mechanism in these sources. 



tron temperature and the spectrum hardens, showing 
signs of becoming saturated at high L. 



4.2. Variable injected power 

Emission at the synchrotron self-absorption fre- 
quency is dominated by the thermal part of the elec- 
tron distribution in case of low injection rate/optical 
thickness (solid line in Fig. 2), but it is dominated 
by the power-law tail in the cases of higher luminos- 
ity. As the injection rate increases, so does the num- 
ber of seed photons for upscattering and the equilib- 
rium electron temperature drops. The Compton y- 
parameter also increases despite the decreasing elec- 



4.3. Variable ratio between magnetic 
and injected compactnesses 

Increasing magnetic compactness increases the elec- 
tron cooling rate, therefore the normalization of the 
cooling power-law decreases, while the thermal dis- 
tribution persists to higher energies (Fig. 3). The 
synchrotron radiation near the self-absorption fre- 
quency is produced by the power-law tail of the elec- 
tron distribution in all three cases. The photon en- 
ergy index a of the thermal Comptonization spec- 
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Fig. 4. Equilibrium electron (left) and photon (right) distributions for different ratios between the power from external soft 
photons (disk) and the injected power: / = (long dashed), 0.1 (short dashed), 1 (dotted), 10 (solid). Other parameters: 
r c = 3.5, -pp = 2, L = 10 37 erg s —1 , r)B = 5. The transition from / = to / = 10 affects the observed spectrum and strongly 
resembles the spectral state transitions in GBH. 



trum f^oci? varies between 0.7 and 0.8. 

4.4. Variable ratio between disk and 
injected power 

Increasing the ratio between the soft disk luminos- 
ity and the injected power leads to the evolution 
from almost purely thermal to a hybrid (mostly non- 
thermal) electron distribution (Fig. 4). The result- 
ing photon distribution changes from hard thermal 
Compton dominated spectrum to the one dominated 
by the disk blackbody, with a weak non-thermal tail. 
This is similar to what is observed in the spectral 
state transitions in GBH [6]. 

5. Conclusions 

1. We have shown that hard electron injection is un- 
able to produce hard Comptonized spectra because 
the resulting high synchrotron luminosity provides 
too many soft seed photons for Compton cooling of 
the electrons. We are therefore able to constrain the 
electron injection mechanism operating in the hard 
state of GBH. 

2. A source at L = 10 37 erg s — 1 , where only syn- 
chrotron and Compton processes operate, produces 
an X-ray spectrum with a « 0.7, which strongly re- 
sembles the hard state of GBH. The feedback from 
the disk does not seem to be needed to produce such 
spectra (in contrast to [7-9]). 

3. At high luminosities L ~ 10 38 erg s^ 1 , in the ab- 



sence of disk radiation, the spectrum is close to sat- 
urated Comptonization, peaking at ~ 5 keV. This 
Wien-type spectrum might be associated with the 
thermally-looking emission in the very high states of 
e.g. GRS 1915+105. 

4. A behaviour similar to what is observed in the 
spectral state transitions in GBH can be reproduced 
by varying the ratio of injected soft luminosity and 
the power dissipated in the corona. 
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